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Abstract 

Schrodinger equation with potentials of the Kratzer plus polynomial type (say, quartic V(r) — 
Ar A + B r 3 + C r 2 + D r + F/r + G/r 2 etc) is considered and a new method of exact construction 
of some of its bound states is presented. Our approach is made feasible via a combination of the 
traditional use of the infinite series tp(r) (terminated rigorously after N + 1 terms at certain specific 
couplings and energies) with several new ideas. We proceed in two steps. Firstly, in the strong- 
coupling regime with G — > oo, we find the exact, complete and compact unperturbed solution of 
our N + 2 coupled and nonlinear algebraic conditions of the termination. Secondly, we adapt the 
current Rayleigh-Schrodinger perturbation theory to our nonlinear equations and define the general 
G < oo bound states via an innovated, triple perturbation series. In its tests we show how all the 
corrections appear in integer arithmetics and remain, therefore, exact. 



I. Introduction 

Bound states in the Coulomb or, after a slight generalization, Kratzer's [Q potential 

V {0} (r) = F/r + G/r 2 

played an important role in the history of quantum mechanics as well as in the various methodical 
studies of its limitations (^] . In applications the potential offered one of the most important exactly 
solvable models of atomic and molecular physics and quantum chemistry. The model is extremely 
transparent, its spectrum of energies E is numbered by integers N = 0, 1, . . . and its wave functions 
are proportional to Laguerre polynomials of degree N @. Almost 30 years ago A. Hautot noticed 
Q that the exact solvability of the Kratzer's potential survives, in a way, its immersion in the 
shifted harmonic confining well, 

V {1} (r) = Cr 2 + Dr + F/r + G/r 2 . 

For any integer N — 0, 1, . . . one may construct N + 1 different "Sturmian" exact wave functions 
proportional to polynomials of degree N. In contrast to the Kratzer's model the Hautot's elementary 
bound states do not form a complete set. In literature the solutions of this type are called quasi- 
exact |J. 

The two classes of forces V{oy(r) and V{iy(r) are usually considered exceptional. This attitude 
finds a deeper mathematical foundation in the representation theory of Lie algebras 0. Still, 
certain remnants of their elementary solvability may be detected in all their further polynomial 
extensions. This has been first noticed by Magyari |8) who generalized the Hautot's method to all 
the polynomial forces 

V {m} (r) = Ar 2m + B r 2 "^ 1 + . . . + F/r + G/r 2 . (1) 

The generalized construction with degree m > 2 did not inspire, unfortunately, too many appli- 
cations. Basically, this was due to its numerous practical shortcomings. In introduction, let us 
only mention that Magyari replaced the linear Schrodinger equation by a nonlinear algebraic sys- 
tem. The solution of these equations seems only feasible for the first two choices of indices m = 



2 



and m = 1, with nonlinearity more or less kept under control. In the present paper we intend to 
demonstrate that such a constraint and a priori scepticism about the feasibility of solution of the 
Magyari's equations are undeserved, for the important Q choice of m — 2 at least. 

The paper is inspired by certain unpublished numerical experiments with Magyari equations 
([[Tol, cf. also Section II). In the first step of a new development we show in Section III that certain 
unsolvable interactions (Q) with m = 2 and with a strong repulsive core G — > oo leave the Magyari 's 
equations exactly solvable at any dimension N. This is our key result. 

In the equally important second step we demonstrate how one can deal with the corrections 
for finite G < oo. Section IV offers a new perturbation method of construction of their sequence 
in closed form. The method emphasizes several close parallels between our nonlinear algebraic 
problem and its simpler linear analogues. Technical details of our Rayleigh-Schrodinger-like new 
perturbation theory are illustrated on our quartic example V{ 2 } (f) in full detail. Section V adds a 
few concluding remarks. 



II. Matrix form of the Schrodinger equation 

An overall methodical background of the Magyari's treatment of the ordinary radial Schrodinger 
equation 



d 2 £{£+!) 



V(r) 



ip(r) = Eijj(r) (2) 



dr 2 r 2 

with angular momentum £ = 0, 1, . . . lies somewhere in between the variational and non- variational 
approaches. While one usually expands wave functions in an appropriate orthonormalized basis, 
Magyari uses the mere Taylor series. His recipe weakens the current emphasis upon methodical 
universality and numerical efficiency. It favors, instead, an enhanced simplicity, analyticity and 
closed, non-numerical form of his quasi-exact bound states in a way which is re-gaining a new 
credit in recent literature 

At the very beginning of our considerations let us note that all the asymptotically polynomial 
forces (pi) with m > 1 imply the same asymptotically exponential form of their bound-state wave 



functions, 

V-(r) = exp {-VA [r m+1 /(m + 1) +Br m /2A + O (r'"" 1 )] } . 

This is an important information. Once we restrict our attention to the mere quartic m = 2 
potential for the sake of definitness and brevity, 

V {2} (r) = Ar 4 + Br 3 + Cr 2 + Dr + Fr- 1 +Gr~ 2 (3) 

we may try to write down the polynomial ansatz 

/ i i \ N 

^ {2} (r) = e X p ar 3 --/3r 2 - 7 rj ^^r"+' +1 . (4) 

At an arbitrary finite integer N > the requirement of its asymptotically correct behaviour fixes 
the parameters a = \[A > 0, = B/2a and 7 = (C — (3 2 )/2a. The next order of asympotic 
analysis also reveals the uniqueness of the N— dependent "termination-admitting" value of the 
coupling constant D = D(N) = —2a(N + I + 2) + 2(3j. This parallels the Hautot's observations 
made at to = 1 H . 

Near the origin the ansatz (^) combines the angular momentum I = 0, 1, . . . with the centrifugal- 
like coupling G > —1/4 into a new quantity I = 1(G) such that G + 1(1 + 1) = 1(1 + 1). After the 
removal of a sign ambiguity in I = -1/2 + ^G + (1+1/2) 2 > -1/2 we may insert our ansatz in 
eq. (^) with potential (^). This transforms the differential Schrodinger equation into equivalent 
recurrences for Taylor coefficients Wj, 

RkUJk-2 + TkUJk-l + SkLOk + PkLOk+l = (5) 

with k — 0,1, . . . , N + 1 and coefficients 

Rn = 2a(N + 2-n), T m = E + 7 2 - /3(2m + 21 + 1) 
(abbreviated as T m = 2T — 2(m —1)0) and 

S m = -2y(m + 1 + 1)- F, P m = ( m + l)( m + 21 + 2) 

(with abbreviations S m == 2S — 2ym and I + 1 = Q). These equations form an over-determined 
linear system of equations for parameters Uj with asymmetric and non-square four-diagonal matrix 
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N + 2 by N+l, 
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Within the approach based on the latter equations one only works with a few exceptional bound 
states. They are elementary and preserve a strong formal similarity to the popular harmonic 
oscillators. Unfortunately, this parallel seems to weaken with the growth of N since the explicit 
solution of our non-linear problem (^) quickly becomes more and more complicated. 

We intend to simplify the equations in question at a cost of their perturbative re-arrangement. 
Such an attempt is motivated by the possibility of a suppression of the most quickly growing 
(binomial) matrix elements at large Q = 1(G) + 1 ~ G 1 / 2 ^> 1, i.e., under the presence of a strong 
repulsive core in potential (|^) . A routine re-scaling of coordinates r — > /1 r in eq. (^) and of the 
related coefficients ui n — > u n = Lu n ^ n in eq. leads to a modified and more transparent matrix 
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form Q u = of the Magyari's eq. (|5|) 



S Q/fi 
aNfi 2 (T S- 2 7 



V 



= o. 



(6) 



\ UN J 

<ur (T-N/3)fi 

In the repulsive fl ^> N regime we notice that |5| 3> N\j\ and \T\ ^> N The two main diagonals 
of Q become approximately constant. 

With the energy E proportional to the new parameter T — T(E) and with the Coulomb coupling 
F contained in S = S(F) a "symmetry" of Q may be enhanced by the choice of the value of the 
scaling parameter /i in such a way that the magnitude of its uppermost and lowest diagonals is 
balanced. For example, with a// 2 = Q/ = t we fix the unique values of /i = (fi/a) 1 / 3 and 
r = (n 2 a)^ 3 . Then we may abbreviate s = S/t and t = fJ.T/r and pre-multiply and split our 
pseudo-Hamiltonian into a sum of matrices Q — + XQ^. Its first component 



Q(°) = Q(°)(M) 



s 1 

t s 
N t 



\ 



V 



t s N 
2 t s 
1 t 



J 



will play the role of our unperturbed Magyari-Schrodinger pseudo-Hamiltonian. In the strong-core 
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limit G — + oo our zero-order Magyari equation reads 



I (0\ 



Q {0) (s,t) 



,(0) 



,(0) 



= 0. 



(7) 



It has the form of an over-complete system of N + 2 equations which should determine all the N + 1 
Taylor coefficients plus the two spectral-like parameters s = and t = f(°). In the purely 
mathematical setting, the difficulties represented by non-square form of their pseudo-Hamiltonian 
Q^(s,t) are quite new and rarely encountered in the mainly linear formalism of the traditional 
quantum mechanics. Of course, the solvability of eq. (0) with m > 2 loses a natural Lie-algebraic 
background and interpretation of its quasi-exact predecessors with m < 1 |tJ . 



III. Unperturbed solutions 

At an arbitrary N, let us introduce the two non-square quasi-unit matrices J and /C such that 
Q(°'(s, t) — Q( o '(0, 0) + sj + tK. Denoting their transposition by a superscript T it is easy to 
imagine that the products J T ' Q^°\s,t) and K7 (s, t) are square matrices. In such a notation 
the necessary and sufficient condition of the nontrivial solvability of eq. (Q) may be formulated as 
a simultaneous disappearance of the two independent secular determinants, say, 



det 



J T Q (0) (0, t) + sl] =0, det \lC T Q (0) (s, 0) + 1 1] =0. (8) 



Our two free unperturbed parameters s — gC' and t = t^°' will be fixed and determined by these 
two coupled polynomial equations at any N in principle. 

Due to symmetries of our generalized, coupled eigenvalue problem (^) the pairs of its (in general, 
complex) solutions will always remain complex conjugate, s^ - 1 = [t^l • For physical reasons we 
must ask whether at least some of these solutions remain real, — . The answer may be 
found by their explicit evaluation. 

At N — the solution of eq. (||) is unique. From the natural normalization u^' — 1 we only get 
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the trivial s = t = 0. At A = 1 the same normalization implies that Uq = — 1/s (e.g., via the first 
row of eq. (0)). The second row defines t = s 2 and our set degenerates to the single cubic equation 
s 3 = 1 . Out of its three different complex roots only one is real and we have the unique physical 
solution s = t = 1. At the next dimension A — 2 an alternative, intermediate normalization 
— 1 helps us to eliminate = —1/s and = — 1/f via the respective first and last row of 
cq. (Q) . This gives s 3 = t 3 and linear relation t = as with the three possible complex constants a 
such that a 3 = 1. One ends up with a triplet of alternative quadratic equations for the unknown s. 
Similarly, one proceeds at the higher integers A. 

The resulting sets of the roots s are sampled in Table |l|. The general pattern of their N— dependence 
is obvious. Once we pay attention to the real solutions only, we arrive at the general formula 



* = 4°+l]( iV ) = *[n+l]( iV ) =iV - 3n ' 



AO) 



0,1,..., [A/2]. 



(9) 



We see that in the limit G — > oo our double eigenvalue problem (g) is solvable in closed form at 
any A = 0, 1,2, . . .. 



A. Pascal-like triangle for the Taylor coefficients 

Taylor coefficients fc(^0i = Q, 1, . .., A define the exact wave functions ip^(r) (||) at each 

physical root (^|). Via a suitable A— dependent normalization one may calculate these coefficients 
in integer arithmetics (i.e., exactly). This is illustrated in Table ^. At the maximal real roots 
s = S[i](A) = N we have 



u^JN) = (-l) w +« 



\ n J 



and our recurrences degenerate to the Pascal triangle for binomial coefficients, 









+ 


fy-A 


[ n J 




{ " ) 




{n-l ) 



(10) 

1 ) 

This implies the elementary form ip(r) = if>m(r) ~ (1 — r//i) N of the related exact wave functions. 
One detects the presence of a degenerate nodal zero of multiplicity A at r = r z = fx. A transition to 
the next, smaller real root s = spi (A) = N — 3 makes the multiplicity of the node in ip( r ) = V'p] ( r ) 



lowered by two. The phenomenon survives iterations in the bracketed subscripts. Our closed formula 
@ for the roots may be complemented by the similar rule N — 2n for the nodal multiplicities. With 
the odd integer N = 2 J + 1 the iterations end at the real spectral root su+i\(2J + 1) = 1 — J. A 
simple, non-degenerate zero appears in the related wave function Wj+i] ( r )- This characterizes the 
first excitation. At even TV = 2K the final choice of the minimal S[x+i\(2K) = —K produces the 
nodeless wave function ipjj^+i]( r )- It describes the ground state. 

A closer inspection of the numerical values of the Taylor coefficients reveals the presence of 
certain non-binomial cases. Most quickly this "anomaly" is spotted in the ground-state coefficients 
fc(2-^0) k — 0,1, ... , 2K. The same pattern re-appears in the first excitations Wrj^-n j(2J + 
1), j = 0, 1, . . . , 2J+ 1 at the higher J for odd N = 2 J + I, etc. Thus, the ground-state set is most 
fundamental and a few more values k (2K) are displayed in Table ||. It definitely confirms a 

surprise. Each of the coefficients proves to be a sum of its three closest upper neighbors, 

Ak + i]A 2K ^ = u m,k( 2K " 2) + u^^K 2) + u$ ]<k _ 2 (2K - 2). 



In an unexpected parallel of the above two-term rule (|10|) our Table |3| forms a new Pascal-like 
triangle. For ground states with different K we have xp^?) (r) = I, ip^) ( r ) = 1 + r I M + 7,2 / M 2 ; 
V , [3?( r ) = (1 + r / fi + r 2 / fi 2 ) 2 and so on. We may generalize this result also to all the wave functions 
V'fe+i] i r ) which correspond to the same (even or odd) integer N = 0,1, . . ., 

<ix](r) = (1 - r/^ N - 2n (1 + r/n + r 2 / M 2 )", n — 0,1, . . . 

This is one of our main results valid, by induction, for all the bracketed subscripts. 



(11) 



B. Roots of the coupled secular equations 

The multiple, seemingly degenerate nodal zeros of solutions (|ll| ) require a "magnification" of their 
vicinity by an appropriate inclusion of a few higher-order corrections 0(l/fl), 0(l/f2 2 ), .... The 
only exceptions are the particular minimal-root ground states 

V$+i]( r ) = (1 + r /» + r 2 /n 2 ) K i N = 2K, s = t = -K, 



!) 



and the first excitations 

i>[j +1] (r) = (l-r/n)(l + r/v + r 2 /Li 2 y, N = 2J + 1, s = t = l-J 

without multiple zeros. In the strong-core phcnomenological regime they may prove useful for 
immediate applications in principle. In such a context even the most simplified schematic example 
V(r) = A r 4 + G/r 2 with G>4 indicates that the forces with large G>1 need not necessarily lose 
a reasonable physical interpretation. A minimum of this function lies at a point tq = (G/2A) 1 / e 
which only slowly moves with G. Whenever necessary, this motion may further be slowed down by 
a simultaneous increase of the coupling A. 

In order to get a feeling for the subtleties of structure of the systems with large G let us pick 
up our ground-state solution with N = 4 and s = t = —K = —2 in s— wave (£ = 0). For simplicity 
let us assume the absence of the Coulombic and cubic terms in our potential V\2}{f)- This means 
F = B = and fixes the values of (3 = and 7 = C/2a with C = 4(A 2 /0) 1 / 3 . At the ground-state 
energy E = — 7 2 — 4(A57) 1 / 3 our special example 

V {2} {r) =a 2 r 4 + 4{a 4 /n) 1/3 r 2 -2a(n + 5)r + n(n-l)/r 2 (12) 

is still only solvable in the limit Q — > 00. Nevertheless, due to our explicit knowledge of the wave 
function, we may now invert the procedure easily. All the higher-order deviations from solvability 
may already be added to eq. (|l^) in an explicit form produced, say, via the insertion of our 
elementary ground state formula for 4'[2]( r ) m the original, differential Schrodinger eq. @). 

For definitness of our illustration let us choose the unit spring-constant coefficient C = 1 at r 2 
in V{2}(r). We get, as a consequence, a strong quartic confinement a 2 — A = \/Ci/8 combined with 
the repulsive linear force possessing an even stronger, 0(£! 5 / 4 ) coupling D — — f2 1 / 4 (f2 + 5)/V2. 
We may re-scale the coordinate r — > rVL 1 ^ in our radial differential Schrodinger equation (|[) with 
force (p^). In the leading-order 1 + 0(1/ ft) approximation this leads to the zero-energy problem 
with the dominant three-term interaction V^ 2 } C 7 ") — const r 4 — const' r + const" jr 2 . All its coupling 



constants are of the same order of magnitude 0(1). It is known |12 that the solutions of such a 



problem may also be elementary and proportional to Laguerre polynomials. 
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We may conclude that in the limit f2 — > oo the whole picture and scheme of elementary solvability 
is nicely self-consistent. We also see the usefulness of our present generalized construction. It 
leads to a richer class of polynomial solutions tp(r) in a fairly nontrivial though still feasible and 
manageable way. What remains for us is to extend the same pattern of self-consistency to all the 
higher-order (i.e., presumably, 1/Sl fe ) corrections. 



IV. Corrections 



At the finite values of the coupling G the strength of the perturbation A = Q — Q^ is measured 
by the parameter A = 1/(/it) = 1/fi. Its smallness is strictly equivalent to our above assumption 
O N. In a perturbative treatment of our bound-state problem Q u — we have to avoid the 
following two most serious obstacles. 



Firstly, our pseudo-Hamiltonian Q(s, t) of eq. 

as well as their three-diagonal difference 

/ 



-7M 1 
-/V -2 7 M 
-2/V 



and its unperturbed simplification Q(°\s,t) 
\ 



3 

-37/1 
-3/V 



N(N-l) 
2 



V 



' ■ • -N'y/j, 
-Nj3^ 2 

are non-square matrices. This means non-linearity, not tractable by the standard perturbation 
formalisms. Secondly, our zero-order pseudo-Hamiltonian itself is a four-diagonal matrix. This 
would cause difficulties even in the linear case, with an unclear idea how one could construct an 
unperturbed propagator. Fortunately, the apparently unavoidable methodical pessimism is not 
in place. With certain care one can proceed in an almost complete analogy with the traditional 
Rayleigh-Schrodinger textbook prescription, 



choosing an arbitrary integer N > 0, 
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picking up a real root = t^ as given by eq. 



postulating, in nonstandard manner, the two different expansions 



= s (o) +As ( 1 )+AV 2 > + 



t = t^ + xt^ + \ 2 t^ + 



plus, more traditionally, 
/ \ 

U 



( (0)\ 



,(0) 



/ ..(X) \ 



+ A 



"0 
,(X) 



,(1) 



(2) 
"0 

,( 2 ) 



.(2) 



yu N j yu N j y u N j y u n J 

• combining all these expansions with input Q = + \ Q^\ 

As a net result we get the hierarchy of equations for corrections. On all the subsequent levels of 
precision 0(X k ), i.e., in the k— th perturbation order with k = 1, 2, . . ., these equations will play the 
role of implicit definitions of the "charge" , "energy" and "wave function" corrections S W, and 
u^ k \ respectively. The set will be initiated by the above zero-order problem (Q). All the subsequent 
new equations will have the same non-square and non-homogeneous common matrix form 



)(0) 



(0, 0) + s (0) J + i (0) K 



7« 



Q«(0,0) + s «.7 + i (1) /C 



+ (V 2 ) J + AC) u (k - 2) + ...+ [s^ J + i« K) u (q) = 0. 



(13) 



Mutatis mutandis, we may parallel the standard textbook perturbation theory of the Rayleigh- 
Schrodinger type [O . 



A. Perturbations of the charges and energies 

Our first important observation is that eq. (^) and properties of our zero-order pseudo-Hamiltonians 
Q(0)( s (0) )t (0)) = Q(o)(o,0) + J + JC imply that these non-square and (N + 2) x (N + 
1)— dimensional matrices always possess the pairs of independent left eigenvectors. Their explicit 
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sample [i?W] T and [-yl !] 7, is displayed in the first part of Table |I| We normalized these vectors of 
dimension N + 2 in the symmetric and antisymmetric manner, respectively. Wherever necessary, 
the Dirac's symbols will be used to denote the similar "longer" columns as "kets" (e.g., = \v^)) 
and the rows as "bras" (say, [z/[ s ]] T = (v^\ etc). 

Marginally, let us admit that there exist alternative possibilities of normalization. The require- 
ment of disappearance of the first or last matrix element in a suitable linear superposition of (v^ s /°^ \ 
may be used to define the "compactified" left eigenvectors (wj = S T J T and (w*\ = ff F IC T '. The 
study of their properties discouraged us from their use in computations since the new doublets do 
not always remain linearly independent. This is exemplified at N — 2 and = 2 in the second 
part of Table |. 

The left action of a left auxiliary eigenbra (ij[ s / J | completely eliminates the vector from eq. 
(ph. We abbreviate 



- (V 2 ) J + K) - ... - (V*- 1 ) J + tc) u {1) 

and arrive at the two independent equations 

(«M | (Ju^)) s<*> + («M | (ICu^ )) t< fc > = («W | S (*-i) ) , 

They may be interpreted as a two-by-two matrix inversion. The symmetric and antisymmet- 
ric normalization of our auxiliary eigenvectors implies that (v^ \{Ju^)) = (v^\(K.u^)) and 
(t>[ a l|( l /u' ')) = — (v^\(K.u^)} . This enables us to invert the left-hand-side matrix in an explicit 



1 «<*> \ ! ( l/(vW\(Ju(°y)) l/(v^\(Ju^)) \ ( (wWjsC*- 1 )) ^ 





I 

1 




~~ 2 

V 



/ 



l/{v^\{JCu^)) 1/( V [ q ]|(/Cm(°))) J \ (v^\E (k -^) 

This is our first final closed formula for corrections s^ k > and v- k > , an extended non-linear parallel to 
the usual Rayleigh-Schrodinger definition of energies. 
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B. Perturbations of the wave functions 



The column vector of corrections is characterized by its re-normalization ambiguity — > 
+ const x u(°\ This follows from the very definition ( |l3| ) and implies that, say, the first 
component of may be chosen as vanishing. Let us indicate such an option by a superscript * 
in — u*( k *> with Uq = for k > 0. In our definition ( |l3| ) this choice of normalization makes 
the pseudo-Hamiltonian matrix (s^°\t^) 7 in effect, lower triangular. In the same equation all 
the known terms may be collected in the single ket 

| r (fc-i)> = | S (fe-D) _ s (fe) jtfo) _ t (k) ^(o) 

entering the right-hadn side of our set of relations (s^-°\ ft ') u*^ — |T^ fe_1 ^). Two of its rows 
(or, more precisely, the two linear combinations of all these rows) have already been used for the 
determination of the quantities and . We are left with the N independent equations. In the 
second decisive step of our method we recommend their choice which drops the last two lines and 
offers the relations 



/ 



R* 



K*0 



\ 



( *(fc-l) \ 



(r(fe-l) 



r N-l 



with 



R* 



1 




s (0) 


2 


t (0) 


s (0) 


N - 1 


t (0) 



V 



<(0) s (0) N 



i.e., an N— dimensional linear matrix equation which defines all the remaining unknown quantities. 
Computationally, it just offers the highly economical evaluation of the higher-order corrections 
via the downwards-running and finite four-term recurrences. 
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An important merit of our proposal is that the main diagonal of our lower triangular (and finite- 
dimensional) matrix R* is safely non-zero. Hence, the matrix itself is regular, deti?* = AH 7^ 0. 
Its lower triangular inversion (-R*) 1 exists and may be interpreted as the present analogue of the 
current Rayleigh-Schrodinger unperturbed propagator. 

The direction of our recurrences is an inessential consequence of the choice of normalization. 
After its alternative specification = u{ k ^ with (u^)n = in all the orders k > we annihilate 
the last column of Q(°)(s(°', t^). Naturally, one omits now the first two rows. Thus, we cut the 
A— dimensional and upper triangular submatrix R± out of Q^°\s^°\ t^) and define the alternative 
unperturbed propagator by inversion, 



V 



,(fc) 



,(fe) 



,(fc) 

l *N 



\ 



( (k-1) \ 



J 



\ T *N-1 J 



Computationally, of course, the higher-order corrections u{ k ^ are again much more easily evaluated 
directly, via the upwards-running four-term recurrences. 



V. Conclusions 

The core of our present paper may be seen in the progress we achieved in the non-numerical 
Magyari-like constructions. By means of their decisive perturbative simplification we were able to 
offer new solutions of certain "next-to-solvable" models. 

In the first step we analyzed the limit G — > 00 of quartic forces V{ 2 }(r) and were able to find 
an exact and complete solution of the underlying nonlinear algebraic equations. This demonstrates 
that far beyond the "safe" territory of the exact and quasi-exact solutions one may still get closed 
formulas much more easily than expected before. 

Another pleasant surprise appeared during our study of corrections (9(1/0). For many years 
we felt deterred by the strong A— dependence of their traditional numerical analyses as well as 
by the quick increase of complexity of their explicit description with m > 2 and dimensions A 
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beyond, say, N = of N = 1 In this context, the simple form of our present quasi-linear 

and N— independent perturbative formalism opens a new direction of development towards many 
practical (say, phenomenological) applications yet to be constructed and appreciated. 

We have demonstrated that the traditional Taylor-series approach [|l5| to the differential Schrodinger 
equation may still offer new surprises. Its main merit lies in its ability of combination of the correct 
threshold behaviour with the correct physical asymptotics of ip( r ) f° r a broad class of elementary 
potentials. As a consequence, many simple models in quantum mechanics may efficiently be clarified 
by means of the related "most natural" power series ansatzs: 

• In place of the differential equation one gets a (sometimes more transparent) set of recurrences; 

• Their unique solution may often be written in a closed (often called Hill-determinant |p^| ) 
form at any trial energy E; 

• The truncated power-series wave functions remain fairly precise in a finite (i.e., most relevant) 
interval of coordinates even after a rough energy guess. 

The transparency of the recipe is often marred by its various practical disadvantages. For ex- 
ample, the efficiency of the Hill-determinant algorithms often lags behind the universal variational 
approaches | fl7j . Also the domains of their applicability may be much more constrained, the un- 



pleasant phenomenon which leads to many non-trivial misunderstandings 18 . Difficulties emerge 
even in numerical implementations of these algorithms: In the fixed-precision computer arithmetics 
one has to avoid the potential instability of recurrences etc |l{| . In such a context the most reliable 
remedy lies, obviously, in an exact termination of the infinite series. 

In our paper we emphasized that the disadvantages connected with the latter strategy are often 
overestimated in applications. We succeeded in showing that the merits of using the algebraic 
Magyari-like equations may often be preserved at a reasonable price. Proceeding in a consequently 
constructive way we have shown that their non-linearity has, in a way, a "weak" form which admits 
a certain "quasi- linear" treatment. 

This attitude is of course trivial for the models of the Coulomb, Kratzer or harmonic oscillator 



1G 



type. For all of them the power series method is able to generate all their exact bound-state solutions 
as terminating. Their exact solvability involves all the angular momenta t and physical energies 
E. Similarly, under weaker assumptions, the same polynomial ansatzs also succeed in offering the 
"quasi-exact" || multiplets of bound states of a finite size N. 

The freedom of choosing the multiplicity N has been significantly extended here. Having noticed 
that the latter two "solvable" categories mean mostly the two- and three-term character of the 
underlying recurrences, we paid thorough attention to their very next, four-term generalization. In 
the context of an enormous interest in quartic oscillators |2^] we imagined that one need not always 
perceive the comparatively short four-term recurrences as prohibitively complicated. At the same 
time, for many phenomenological applications, they offer a partially non-numerical insight into the 
new class of interactions which is already fairly rich. 

In the light of our present results, the nonlinearity of the underlying algebra need not represent a 
serious technical obstacle. This is our main message. Via a detailed analysis of our present quartic 
polynomial examples we have shown that the explicit solution of the related Schr"odinger and 
Magyari equations remains feasible in an almost complete parallel with their quadratic predecessors. 
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Table 1: The key to our nonlinear problem (||): Its roots — for the first few N ^= 0. 



N 


1 2 3 


2 Re s(°) 
2 Im s(°> /V3 


2 -1 4 1 -2 -2 6 3 -3 -3 
±1 ±1 ±2 ±1 ±2 ±1 ±3 


iV 


4 


2 Re s(°) 
2 Im s(°) /V3 


8 5 2 2 -1 -1 -4 -4 -4 
0±10±2±1±3 0±2±4 


iV 


5 


2 Re fl (°) 
2 Im s(°) /V3 


10 7 4 4 1 1 -2 -2 -2 -5 -5 -5 
±1 ±2 ±1 ±3 ±2 ±4 ±1 ±3 ±5 
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Table 2: Real eigenvectors of our pseudo-Hamiltonians Q(°\ 



N «(°) 


u 0^ u 2°' u 4°' u 6^ 







1 1 




2 2 
_1 


1-21---- 
1 1 1 


3 3 



1-3 3-1 

10 0-1--- 


4 4 
1 

-2 


1-46-41 

1 -1 -1 1 - - 

1 2 3 2 1 - - 


5 5 
2 
-1 


1 -5 10 -10 5 -1 
1-2 1-1 2-1 
1 1 1-1-1-1 


6 6 
3 

-3 


1 -6 15 -20 15 -6 1 
1 -3 3 -2 3 -3 1 
1 -2 1 
1 3 6 7 6 3 1 
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Table 3: Generalized Pascal triangle (ground states: N = 2K, s\k+i\ = —K). 



K 


(0) 


(0) 
U K-3 


(0) 
U K-2 


(o) 

U K-1 


(o) 


(o) 

U K+1 


(0) 
U K+2 


(o) 

U K+3 


(o) 

U K+4 

















1 














1 











1 


1 


1 











2 








1 


2 


3 


2 


1 








3 





1 


3 


6 


7 


6 


3 


1 





4 


1 


4 


10 


16 


19 


16 


10 


4 


1 


5 


5 


15 


30 


45 


51 


45 


30 


15 


5 


6 


21 


50 


90 


126 


141 


126 


90 


50 


21 
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Table 4: Pairs of independent left eigenvectors of Q^(s, s). 





sym. 


and antisym. 


v k 


shortened 




and 


overlaps 




k = 


o 


1 


2 


3 4 


5 


o 


1 


2 


3 4 




AT 


s 




















\ i \^ / / 


U 


U 


i 


i 








-i 
1 








1 






-i 
i 


-i 








-i 
1 








1 


-i 
1 


1 


1 


o 
-z 


-i 
1 






-i 
1 


-1 






O 






1 





-1 






- 1 


1 






1 
± 


o 
z 


2 


2 


-1 


-1 


z 




U 


1 


-1 




Q 
O 






a 
U 


1 


-1 


A 
U 












o 




-i 

-1 


1 


1 


1 


1 
1 




o 
O 


1 


9 




o 






3 


-1 


1 


o 

-3 




o 

z 


1 
1 






o 


3 


3 


1 


1 


-2 


1 1 




1 





1 

- 1 


1 
1 


Q 

-y 






i 
i 


i 

-i 


A 
U 


i i 
i - 1 




1 


1 

- 1 


n 
u 


1 
1 


o 







2 


-1 





-1 2 




2 





o 


-1 


3 






2 


1 





-1 -2 




- 1 








2 


3 


4 


4 


1 


-2 


1 


1 -2 


1 


1 


-1 





1 -1 


9 






1 





-1 


1 


-1 


- 1 


1 





-1 1 


-9 




1 


7 


1 


-2 


-2 1 


7 


7 


-3 


-1 


-1 4 


9 






1 


-1 





1 


-1 


4 


-1 


-1 


-3 7 


3 




-2 


1 


1 


1 


1 1 


1 


3 


1 


2 


1 2 


9 






3 


-1 


1 


-1 1 


-3 


2 


1 


2 


1 3 


3 
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